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ABSTRACT 

Models invoking large populations of obscured AGN are known to provide good 
fits to the spectrum of the X-ray background and the observed soft and hard X-ray 
number counts. An important consequence of these models is that significant quantities 
of dust are required in order to provide the obscuring medium, which will be heated 
by the nuclear radiation from the AGN, and will radiate in the thermal infra-red. We 
therefore model the properties of the dust, and predict the contribution of obscured 
AGN to the intensity of the far infra-red background and the sub-mm source counts 
and redshift distribution, in order to ascertain whether our models are constrained 
further by the data available at these wavelengths. Our conservative models predict 
a contribution of between 5 and 15 per cent of the far infra-red background intensity, 
whereas if we use more extreme values for several parameters, this value may reach 
33 per cent. This suggests that AGN may only form a significant minority of sub- 
mm sources, in agreement with the results of spectroscopic follow-up. Our models 
thus show that there is no inconsistency between obscured AGN models for the X- 
ray background and either sub-mm source counts or the intensity of the far infra-red 
background. We further propose that obscured AGN may explain the sub-mm emission 
associated with Extremely Red Objects (EROs). Finally, we make predictions for the 
redshift distribution of AGN sources detected in faint sub-mm surveys, which should 
allow future tests of this obscured AGN model. 
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1 INTRODUCTION 

Recently, enormous advances have been made in the field of 
sub-mm astronomy, most notably with the commissioning 
of SCUBA, the Sub-mm Common User Bolomet er Array 



( |Gear fe Cmmiiigliam 1994 | HuUand el al. 1999 ] ), un the 



James Clark Maxwell Tclc3copc (JC MT). The 850/x m sur - 
veys of Smail, Ivis on & Blain (1997), |Hughes et at| (1998) 



and Barger et al. (1998) suggest that the sub-mm counts 



show strong evidence of evolution. The sources detected at 
this wavelength appear to be dusty, high redshift galaxies. 
However, the origin of the emission is not clear, and these 
distant dusty galaxies could be powered by either starburst 
or AGN activity. The question of which component dom- 
inates is highly important: if starbursts dominate, then it 
would imply that a significant amount of star-formation has 
been missed by studies which use ultra-violet techniques 
to claim that the star-formation rate (SFR) decreases in 
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galaxies at high redshift (Madau et al. 1996). Metcalfe et al 
(1996) have previously suggested that models which include 
dust with a monotonically increasing SFR towards higher 
redshift can als o give good fits to faint galaxy counts and 



colours (see also JWang^P991 j_ Sronwall , 
fc Shanks 1997| ; jgtcidcl ct al. 1999[ ). 



Koo 199! 



Campos 



It has also been known for some time that some ultra- 
luminous infra-red galaxies (ULIRGs), which are strongly 
star- forming, dusty galaxies, also contain buried QSO nu- 
clei. Images of ULIRGs in polarized light have shown highly 
anisotropic structure, such as that observe d in the IRAS 
galaxy, F10214+4724, ( [Lawrence et al. 1993 ). Similar polar- 
ization structure is seen in Seyfert 2s and high redshift radio 
galaxies, which is thought to be indicative of non-uniform 
illumination consistent with the Unified Model of AGN. It 
has been inferred that the QSO and ULIRG phenomena 
are closely related, with ULIRGs being postul ated as be- 
ing Type 2 objects, or "QSO-2s" (|rlines 199^ . Genzel et 
al. (1998) have made ISO observations of ULIRGs in order 
to determine which is the dominant emission mechanism at 
FIR wavelengths. They find that ~ 75 per cent are powered 
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by star-formation, while the remaining ~ 25 per cent are 
AGN (see also Rigopoulou, Lawrence & Rowan- Robinson 
1996). A massive dusty torus will also contain huge amounts 
of molecular gas, and therefore by definition will also be a 
prime site for star-formation. It has become increasingly evi- 
dent that both AGN and starburst activity are often present 
in the same object, and therefore it is naive to assume that 
the two processes are either independent or mutually exclu- 
sive. 

The first spectroscopic follow-up of the faint sub-mm 
sources tends to show a mixture of dusty starburst and 
AGN. However, in obscured sources, it is frequently diffi- 
cult spectroscopically to disentangle starburst from AGN 
components. Thus, in a SCUBA survey of four distant clus- 
ters of galaxies (Smail, Ivison & Blain 1997), the brightest 
source, SMM 02399-0136, was fou nd to be a hyperluminous, 
active galaxy at redshift z = 2.8 (Ivison et al. 1995). Frayei 
et al.|(1998) compare the measured Lfi r/I/qq ratio of thi s 
high redshift galaxy, SMM 02399-0136 ( [ivison et al. 1998| ), 
to that of the local starburst, Arp220, and find that it is 
twice as high. Since the FIR emission from Arp220 is purely 
from merger-induced star-formation, they infer that approx- 
imately half of the FIR emission from SMM 02399-0136 must 
therefore be due to a dust-enshrouded AGN. Further exam- 
ples include the ultraluminous BALQSO, APM 08279+5255, 
which at z ~ 3. 9 is apparently th e most luminous object 
currently known ( Lewis et al. 199S ) , and the sub-mm source 
B 1933+503, t hought to be a high redshift (z > 2) dusty 
radio quasar (Chapman et al. 1999). Such surveys suggest 
that starburst galaxies dominate over AGN at faint sub-mm 
fluxes, but the area of sky and number of sources detected 
is still quite small. 

Of course, large populations of such dust-enshrouded 
QSOs and AGN have been hypothesized in order to explain 
the observed shape of the X-ray backgrou nd (XRB) spec- 



tru m (Madau, Ghisellini fc Fabian 1994; Comastri et al 



1995 



Ourm & Shanks 1999). In addition, the hard X-ray 



source counts are considerably highe r than at soft X-ray en- 
ergies ( Georgantopoulos et al. 1997), implying the presence 
of an additional hard-spectrum population to the broad- 
line QSOs known to make up the bulk of the sources at 
lower energies. However, one of the central features of ob- 
scured QSO models is the presence of significant quantities 
of gas and dust surrounding the active nucleus, providing the 
means of absorbing a large fraction of the intrinsic radiation 
and affecting the observed properties in the optical and X- 
ray regimes. But this energy must escape somewhere, and 
through heating of the dust, the absorbed flux is re-radiated 
at far infra-red (FIR) wavelengths, providing an important 
test of the models. The amounts of dust invoked are quite 
considerable, and combined with the huge quantity of energy 
that must be radiated, plus the wide redshift distribution of 
these sources, would imply a substantial contribution to the 
FIR/sub-mm background radiation and source counts. This 
motivates our investigation here into the impact of obscured 
QSO models, used initially to explain observations at X-ray 
energies, at much longer wavelengths. 

While this paper was in preparation, Almaini, Lawrence 
& Boyle (1999) have published sub-mm count predictions for 
obscured QSOs which give similar results to those described 
here. However, the conclusions of these authors are based on 
a simpler model where both the fit of the obscured QSOs to 



the XRB and the obscured QSO sub-mm spectrum have to 
be assumed. We believe that the sub-mm predictions pre- 
sented here are more robust since they are based on mod- 
els which have been shown to give good fits to the X-ray 
background over its full energy range. In addition, in our 
method the absorption of the X-ray and optical spectra to 
produce the sub-mm spectra leaves no question as to their 
self-consistency. 

In Section ^, we describe the latest measurements of the 
spectrum and intensity of the far infra-red background, and 
recent observations of the sub-mm source counts. Section ^ 
discusses how our obscured AGN model can be extended to 
these longer wavelengths, and describes how the emission 
from the obscuring dust torus has been modelled. The re- 
sults are presented in Section [|, and are compared with the 
observed sub-mm number counts and the intensity of the 
far infra-red background at these wavelengths. We discuss 
the implications of our results in Section H, and make pre- 
dictions for the number-redshift relation for future 850/im 
surveys in Section []. 

So me of these results have already been discussed by 



Gunn (1999). 



2 FAR INFRA-RED AND SUB-MM 
OBSERVATIONS 

It is extremely difficult to measure the extragalactic far 
infra-red background (FIRB), due to the presence of fore- 
ground components from interplanetary zodiacal dust emis- 
sion (peaking at ~ 20/xm) and interstellar dust emission 
from our Galaxy (peaking at ~ 150pm). This zodiacal and 
galactic contamination must be carefully modelled, and sub- 
tracted to leave the extragalactic background. In addition, 
the cosmic microwave background (CMB) also contributes 
at these wavelengths, peaking at ~ 1 mm, and this must also 
be accounted for. Once these components have been accu- 
rately modelled and removed, it should then be possible to 
detect any extragalactic FIRB. 



The first detection of the FIRB was claimed by Pugct 



et alj (1996), who used data from FIRAS on board CO BE, 
taking advantage of the FIR window from 200 — 800/^m, 
between the peaks of the zodiacal emission and the CMB, 
and found that the intensity of the extragalactic background 
has the form: 



vl v ~ 3.4 



A 



400/im 



x 10" a W ra- 



in the range 400/im < A < 1000/^m. More recent measure- 
ments have been made at shorter wavelengths using DIRBE, 
also on board COBE (Schlegel, Fink beiner & Davis 1998; 



Hauser et al. 1998; Fixsen et al. 1998) 



At longer wavelengths, there are several surveys cur- 
rently in progress with the aim of resolving the source pop- 
ulations contributing to the sub-mm background. These sur- 
veys take advantage of the increased sensitivity and resolu- 
tion now available with the SCUBA camera on the JCMT. 

SCUBA sub-mm surveys can be divided into tw o types, 
the fi rst be ing pointed observations of blank f ields (Hughet 



st al.| 1998; |Barger et al. 199$ [Eales et al. 1999, ). The second 
strategy is to make pointed observations of clusters of galax- 
ies, in order to take advantage of the gravitational amplifica- 
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tion due to the lensing mass of the cluster (Smail, Ivison & 
Blain 1997) . In this way, sources can be detected that would 
be fainter than the flux limit possible without the amplifica- 
tion factor, and in addition, due to the magnification of the 
source plane, there are fewer problems with source confu- 
sion. The 850/^m number counts currently published can be 
fitted with a power-law of N{> S) = 7.9 x 10 3 S _1 - 1 , where 
N is the number of so urces per square de gree detected above 
a flux limit of SmJy (Smail et al. 1998). 



These recent measurements of both the spectrum of the 
far infra-red background, and the sub-mm number counts, 
now provide strict constraints with which to test theories of 
galaxy evolut ion, star-formation history, etc. In the SCUBA 
Lens Survey ( Smail et al. 1998 ) sample of sub-mm sources, 
at least 20 per cent were fo und to show evidence for AGN 
activity (Barger et al. 1999). In the next Section, we inves- 
tigate the implications of the obscured AGN hypothesis, by 
modelling their properties at sub-mm wavelengths. In Sec- 
tion [|, we predict the number of such sources expected from 
our models compared with the data described above, and 
put limits on their contribution to the far infra-red back- 
ground. 



Flat distribution of colu 
Tilted distribution 




20 



22 24 
log N H (cm ) 
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Figure 1. Comparison of the two distributions of obscured QSO 
populations used in this paper, where the size of each population 
is normalized relative to the unobscured \og(N[i / cm - 2 ) ~ 19.5 
population. Our model of a flat distribution of columns is shown 
as hatched regions, and the tilted distribution is shown as cross- 
hatched regions, where the two models are shown slightly offset 
from one another for clarity. The tilted distribution has a greater 
proportion of high column density objects, as used to improve the 
fits to the data in a go = 0.5 universe (Gunn & Shanks 1999). 



3 MODELLING 

He re we aim to us e the obscured AGN models developed 
in 3unn & Shanks (1999) in order to make self-consistent 



predictions at sub-mm wavelengths. We assume that the lu- 
minosity absorbed at high energies is then reprocessed and 
emitted at lower energies, with a known thermal spectrum. 
The fluxes expected from such objects can then be esti- 
mated, from which we can predict sub-mm source counts 
and their contribution to the far infra-red background. The 
assumptions we have made are the following: 

• The intrinsic 0.3 — 3.5 keV X-ray luminosity of each 
source is known, and is assumed to come from the zero- 
redshift X-ray luminosity function, with parameters as de- 
scribed in Table 

• The column density, Nh, perceived by the optical, 
ultra-violet, and X-radiation is assumed to be a measure of 
the intrinsic amount of dust present, and is not affected sig- 
nificantly by the viewing angle. The covering factor, f C ov, of 
this obscuring material is defined to be the fraction of lines 
of sight for which a constant Nh is seen, and the remaining 
lines of sight are unobscured. 

• All the absorbed flux goes into heating up the dust and 
gas in the obscuring medium, whatever the assumed geom- 
etry, which is then re-radiated isotropically in the thermal 
infra-red. 



3.1 Column density distribution 

Several different column density distributions have been pro- 
posed in order to explain the observed spectrum of the X- 
ray background. For example, Madau et al. (1994) used a 
bimodal model, consisting of a population of unabsorbed 
sources plus 2.5 times as many absorbed sources whose 
column densities had a Gaussian distribution with mean 
log(Af«/cm- 2 ) 



four sub-classes of absorbed populations, each with a differ- 
ent normalization with respect to the unabsorbed popula- 
tion. 

Here we use two models, the first b eing a flat distri- 



bution of columns, which was found in Gunn & Shanks 
(1999) to give good agreement with the observed X-ray 
source counts and XRB spectrum for a qo = 0.0 universe. 
The second model has a tilted distribution, with a greater 
proportion of high columns than in the flat model, which 
gives better agreement with the data in a qo = 0.5 uni- 
verse. The differences in the observed properties of each 
population of sources are attributable solely to the column 
density. The flat distribution uses seven populations evenly 
spaced in log space: \og(N H / cnT 2 ) = 19.5, 20.5,21.5,22.5, 
23.5, 24.5, 25.5, each containing the same number and intrin- 
sic luminosity function of AGN. The tilted distribution uses 
six populations, with the relative normalizations obeying the 
relation: 

$(tfH) = {l+0.5]og(^£)}*'. 

The two distributions are compared in Fig. |lj 



3.2 X-ray luminosity function and evolution 

Since we are assuming that the observed sub-mm emission 
results from reprocessed nuclear X-ray and optical radiation, 
we take as a starting point a two power-law zero-redshift X- 
ray luminosity function (XLF): 



$x(Lj 



^ X ^ A44 



7 2 r *(72-7l) 
. Li \ 



" -^-44 -^-44 



L x <L* x (z = 0) 
L x > L* x (z = 0), 



24.0 ± 0.8. Comastri et al. (1995) used L* x {z) = L* x (0) 10 



where $* is the normalization of the XLF, and Lx 44 is the 
0.3 — 3.5 keV X-ray luminosity in units of 10 44 erg s _1 . Here 
the luminosity evolution is parametrized as either polyno- 
mial evolution: 

(7zZ+7zZ 2 ) 
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Table 1. Parameters of four Boyle et al. (1994) model fits to the ROSAT and Einstein EMSS QSO X-ray luminosity function and its 
evolution, as used in our calculations. The first column contains our adopted descriptive model name, POW(go) or POL(go), with the 
corresponding model name from Boyle et al. (1994) in the column headed B94. The XLF is parametrized by a broken power-law, with 
indices 71 and 72 below and above the break luminosity L* x . The evolution is described by either power-law or polynomial evolution, in 
both qo = 0.0 and go = 0.5 universes. 



Model 


B94 


go 


Evolution 


71 


72 


logZ^(0)t 


7z 


Zcut 


7i 




POW(O.O) 


G 


0.0 


(1 + 


Z)"l* + Z c ut 


1.53 


3.38 


43.70 


3.03 


1.89 




0.79 


POW(0.5) 


11 


0.5 


(1 + 


Z) 1 ' + Z cut 


1.36 


3.37 


43.57 


2.90 


1.73 




1.45 


POL(O.O) 


K 


0.0 




+ <z 2 


1.50 


3.35 


43.71 


1.14 




-0.23 


0.84 


POL(0.5) 


L 


0.5 


~1zZ 


+ <z 2 


1.26 


3.32 


43.52 


1.15 




-0.26 


1.89 



t where L* X (Q) is the 0.3 - 



t where <fr* x is in units of 10 



3.5 keV luminosity in units of ergs . 
~Mpc" 3 (10 44 ergs" 1 )" 1 . 



or power-law evolution: 
L x (z)=L* x (0)(l + zy*, 

where a maximum rcdshift. z„„> . at which the evolution stops 
is incor borated, such that: 



L*x(z) = L x {zcut) 



Z ^> Z C u 



The parameters for the XLF and its evolution have been 
taken from the fits to ROSAT and Einstein EMSS QSO 
number counts by Boyle et al. (1994), and are listed in Table 
|lj The broad line QSOs used to define the XLF correspond 
to our populations with log (JVh/ cm~ 2 ) = 19.5 and 20.5, so 
we normalize the populations accordingly. 

3.3 Canonical X-ray/optical QSO spectrum 

The intrinsic QSO spectrum, F(E), is assumed here to con- 



sist of two power-law s, with spectral in dices a x = 0.9 (Nan- 
dra fc Pounds 1994| ) and a opt = 0.8 ( |Francis 199^ ) at X- 
ray and optical energies respectively. In addition, the X- 
ray power-law is modified by the effects of reflection. The 
relative normalization is defined by a power-law of slope 
ol ox = 1.5 joinin g 2 keV and 2500 A ( Tananbaum et al 



197E; Yuan et al. 1998). The behaviour of the spectr um be- 



tween these two regimes is not well known however. Zheng 
ot al. (1997) investigate the far ultra-violet properties of 
a sample of high redshift QSOs and find that the radio- 
quiet QSO spectrum can be approximated by a broken 
power-law, where the spectral break occurs at ~ 1050 A 
(~ 0.01 keV). The spectral index at longer wavelengths is 
a op t ~ 0.86±0.01 (c.f. a op t ~ 0.8 used her e), steepening sig - 
nificantly to shorter wavelengths (see also Laor et al. 1997). 
For simplicity, we therefore assume that the spectral break 
can be approximated by a step discontinuity at 0.01 keV. 



Nh ~ 10 24 cm -2 due to photo-electric absorption alone 
(see Fig. ^|), we do not take electron scattering into ac- 
count here as it does not affect our res ults significantly. The 



photo-elec tric absorption coefficients (Morrison fc McCam 



mon 1983 ) and can be evaluated using xspec ( /Vma.ud 1996 ) 



In t he optic al, the dust extinction laws from Howarth (1983) 
and Beaton (1979) are used. A constant gas to dust ratio is 
implicitly assumed here. As the column increases, so does 
the amount of the continuum which is destroyed by these 
processes. 

The absorbed luminosity can be approximated by cal- 
culating the energy at which the optical depth is unity for 
each process, E(t x = 1) and E(tau op t = 1) respectively, 
as a function of the column density, and then by assuming 
that all the luminosity emitted between these two energies 
is absorbed. Using the relationship t(E) — g{E)Nh, the 
cross-section scales as a(E) = 1/Nh, for an optical depth of 
unity. The absorbed luminosity is then calculated from: 

L^bs = fcov / F(E) dE. 

J E( Topt = l) 

The range of intrinsic luminosities is defined by the X- 
ray luminosity function (XLF), so the absorbed luminosity 
is normalized by the known X-ray luminosity, Lx- 



Lx = 



F{E)dE. 



The fraction of the luminosity which is absorbed can then be 
calculated. Fig. ^| shows how the absorbed fraction increases 
with column density, saturating above Nh ~ 10 25 cm~ 2 , 
as no radiation whatsoever can escape unaffected from the 
nucleus, and all the emitted energy goes into heating up the 
obscuring medium. 



3.4 Absorbed X-ray/optical QSO spectrum 

The next step is to estimate the absorbed luminosity, L a b s , 
responsible for the heating of the torus, for each column den- 
sity used. At X-ray energies, the opacity is dominated by 
photo-electric absorption at columns of Nh < 10 24 cm -2 . 
Above Nh ~ 10 24 cm~ 2 , the obscuring medium becomes 
Compton thick due to electron scattering, such that the 
effective optical depth is T e g = r ph + r os . Since the ab- 
sorbed fraction of the luminosity is almost unity above 



3.5 Far infra-red QSO spectrum 

Finally, we assume that all the absorbed radiation, L a b s , 
has to escape as thermal emission from the dust in the far 
infra-red, and therefore that Lfir = L a b s . We also assume 
that the FIR luminosity is isotropic, and therefore that the 
received flux is independent of viewing angle. For the case 
for which the obscuring medium is also isotropic, this is a 
realistic assumption, as self-shielding of the inner regions 
means that only radiation from the outer, cooler layers of 
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Figure 2. The fraction of the total X-ray and optical luminosity 
absorbed, as a function of the column density, Nh , for a covering 
factor of unity. Once the column is as high as ~ 10 25 cm -2 , 
effectively all the energy is absorbed. For f cov = 0.5, the fraction 
absorbed will be half that shown here. 



dust will be received. However, for a toroidal geometry (e.g., 
fcov ~ 0.5), radiation from the hot dust in the innermost re- 
gions will be able to escape from the top of the torus, and 
for face-on viewing angles, the spectrum will be broadened 
due to the superposition of components at a range of differ- 
ent dust temperatures ( Pier fc Krolik 1992j ). This has the 
effect of boosting the flux at short wavelengths (< 10/im) 
but does not change the spectrum significantly at the wave- 
lengths at which we are primarily concerned (> 100/im), as 
the hot emission component makes little contribution, and 
the torus is optically thick. As the true geometry is still un- 
clear, neglecting this effect means that the predictions from 
our models at short wavelengths should be taken to be lower 
limits. We therefore choose to approximate our obscuring 
medium by isothermal dust, emitting isotropically. 

By convention, at these energies, frequency units are 
used, and we write: 



P(v e ) dv e 



Assuming optically thin dust emission, the Planck func- 
tion, B(v,T), is modified by an opacity law, where the opac- 
ity depends on both the dust grain composition and the 
size and shape distribution of the gr ains, and which can be 
parametrized as n d oc i/ 13 . Following Cimatti et al. (1997), 
we use the opacity law: 



K d = 0.15 



( v ° v 

1 250 GHz/ 



cm g 



The dust te mperature is take n to be in the range 30 K < 



T d < 70 K ( [Haas et al. 199Sj ; |Benford et al. 1998[ 1999). 
Since the emitted FIR luminosity is constrained by the X-ray 
luminosity, the effect of increasing the assumed temperature 
means that the total normalization must be reduced in order 
to keep I/fir, constant, and vice versa. 

The emitted power, P(v e ), can therefore be 
parametrized as follows: 



P{v e ) = 4nK d (v e )B(v e ,T d )M d , 



(1) 



To calculate the received flux from such a source, we 
use the relationship: 



P(u e ) 

b(v ) du = - — — du e 



(1+2 



dv , 



where Dl is the luminosity distance and dv e 
to give: 



S(u ) 



(1 + z)K d (ve)B(v a , T d )M d 



(1 + z) dv , 



(2) 



To calculate the integrated source counts as a function 



of flux d ensity, we then use the same method as in 



Gunn & 



Shanks (1999), starting from the 0.3 — 3.5 keV XLF, and us- 



ing the above relationships between Lx and S{v ). Note that 
in the sub-mm regime, traditionally flux density is used, with 
units of Janskys, (Uy = 10~ 23 erg cm -2 s _1 Hz -1 ), rather 
than the broad-band flux used at X-ray energies (units: 
ergcm _2 s _1 ), or optical magnitudes. 



4 SUB-MM PREDICTIONS 

Here we investigate the contribution to the sub-mm source 
counts and the far infra-red background predicted by our 
obscured AGN model. We look at the effects of changing 
certain parameters, such as the covering factor and temper- 
ature of the obscuring medium, the luminosity evolution of 
the AGN, and the maximum redshift at which these sources 
exist. At present, these parameters are not well constrained, 
particularly at high redshift. However, despite eviden ce that 
the space density of QSOs declines beyond z ~ 3 (Shaver 
et al. 1996), we know that QSOs exist at high redshift. New 
search techniques are discovering more such objects all the 
time, for instance the three QSOs found recently by the 
SDSS Collaboration (Fan et al. 1998), all with redshifts in 
the range 4.75 < z < 5.0. As we are using luminosity func- 
tions and evolutionary models determined from X-ray se- 
lected QSOs, these have diverging properties above z ~ 2, 
and these can therefore be taken to span the range of likely 
properties. By invoking the most extreme cases, we are able 
to put firm upper limits on the obscured AGN contribution 
to both the source counts and FIRB intensity. 



4.1 The covering factor of the absorbing material 

In order to remain c onsistent with the assumed flat column 
distribution used in Gunn & Shanks (1999), we must also 



consider the covering factor of the obscuring material, as this 
affects the intrinsic distribution of column densities. Here, 
we take two examples, first that the obscuring material is 
isotropic, i.e., with a covering factor f cov = 1, and secondly 
that a torus covers half the sky as perceived by the nucleus 
for all sources, i.e., f C ov = 0.5. We assume that these two 
values will span the true range of covering factors. Two in- 
trinsic distributions giving rise to a perceived flat column 
distribution, as used in the qO = 0.5 model, are shown in 
Fig. | 

A method for differentiating between these two pro- 
posed geometries for the obscuring material in AGN, is to 
determine the fraction which are highly luminous at sub- 
mm wavelengths. By definition, very little dust exists along 
the line-of sight to broad-line QSOs. Therefore, if the ab- 
sorbing medium is isotropic, the dust content of broad-line 
QSOs must be intrinsically low, with low sub-mm emission. 
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Table 2. The contribution to the intensity of the FIRB at 850/^m predicted by our obscured AGN model. The parameters used are as 
described in the text, where we have investigated the effects of changing the X-ray luminosity function and evolution, the covering factor 
of the obscuring medium f cov and its temperatu re T^, and the maximum redshift z max - The predictions are compared with the observed 
850/xm background intensity from Fixsen et al.| (1998) of /fire = 5-03 X 10~ 10 Wm~ 2 sr _1 . 
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Figure 3. Comparisons of two different intrinsic column density 
distributions of obscured QSO populations, which give rise to the 
same perceived flat distribution, as used in the qO = 0.5 model. For 
an isotropic absorber, the covering factor is unity, i.e., f CO v = 1, 
and therefore an intrinsic flat distribution will also be perceived as 
flat (hatched regions); however, for an absorber for which f cov = 
0.5, a larger number of obscured sources are required, as half 
will be observed to be unobscured (cross-hatched regions); the 
two distributions arc shown slightly offset from one another for 
clarity. For the tilted distribution of columns used in go = 0.5 
scenarios, a similar correction must be made when f cov = 0.5 is 
assumed. 



However, f or a toroidal str ucture in the spirit of the Uni- 
fied Model (Antonucci 1993), broad- line QSOs could contain 



large amounts of dust in a plane perpendicular to the line 
of sight. This is consistent with the large infra-red bump 
observed in the 3000 A to 300/im spectru m of a sample of 
Palomar Green QSOs (ganders et al. 198S| ). Therefore, if all 
QSOs contain large quantities of dust, they will be strong 
sub-mm sources. Proposed sub-mm observations of X-ray 



selected QSOs using SCUBA/JCMT will shed further light 
on this question in the near future. 



4.2 The effects of changing the XLF parameters 

We first investigate the predicted source counts for the ob- 
scured QSO model with a fiat distribution of columns, inte- 
grated over redshifts < z < 5. The X-ray luminosity func- 
tion used is evolve d according to the power-law and polyno- 
mial models from Boyle et al. (1994), for both go = 0.0 and 
qo = 0.5 cosmologies. The parameters of these models are 
detailed in Table hi 

First, we take as our fiducial model a power-law pre- 
scription for the luminosity evolution, and a qo — 0.0 cos- 
mology, which we denote as POW(O.O). This model will be 
used hereafter, unless stated otherwise. In Fig. ^, we show 
the source counts predicted by this model, compared with 
the SCUBA data described in Section | and Fig. |. A cov- 
ering factor of f cov = 1 is assumed, for a dust temperature 
of Td = 30 K, integrated over < z < 5. The contribu- 
tion from each population of obscured QSOs is shown by a 
dashed line, with the total source counts denoted by a solid 
line. The reverse trend to that observed at X-ray energies 
is seen here: the populations with the smallest column den- 
sities give the lowest sub-mm fluxes, with the high column 
density populations dominating the sub-mm source counts 
due to the large dust masses present. This model predicts 
~ 20 per cent of the source counts at 2 mjy, and flattens off 
at fainter fluxes, but will provide ~ 17 per cent of the FIRB 
at 850/im. 

In Fig. ||, the predictions for all four Boyle et al. mod- 
els are shown, for covering factors of (a) f cov = 1 and 
(b) fcov = 0.5, using the intrinsic column distribution that 
each covering factor implies. As expected, the low qo mod- 
els give a much larger contribution to the source counts 
than the go = 0.5 models, due to volume effects. Tak- 
ing the measured value for the 850/an FIRB intensity of 
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Figure 5. Total predicted 850/rni source counts of obscured QSOs for models using power-law and polynomial evolution, go = 0.0 and 
go = 0.5 cosmologies, compared with the observed sub-mm counts as described in Fig. ^. A covering factor of (a) f cov = 1, and (b) 
fcov = 0.5, is assumed, with a dust temperature of Td = 30 K, integrated over < z < 5. These models can account for between 17 per 
cent (POW(O.O), fcov = 1) and 4 per cent (POL(0.5), fcov = 0.5) of the FIRB. Note the change of scale with respect to Fig. 0. 
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Figure 4. Predicted 850/^m source counts of obscured QSOs, 
compared with the observed counts as described in Section 0. The 
filled circles show the results of the SCUBA Lens Survey (Blain 
et al. 1999), and the other points are as labelled: S97 - Smail, 
Ivison & Blain (1997); B98 - Barger et al. (1998); H98 - Holland 
et al. (1998); E99 - Eales et al. (1999); HDF, P(D) - Hughes 
et al. (1998). A dust temperature of Td = 30 K and a covering 
factor fcov = 1 is assumed, integrated over < z < 5, for power- 
law luminosity evolution and go = 0.0, model POW(O.O). The 
contribution from each individual population of obscured QSOs is 
shown as a separate dashed line, marked with the column density, 
and the total is shown by the solid line. In contrast to the X-ray 
number counts, here the biggest contribution comes from the most 
highly obscured objects, since these contain the largest amounts 
of dust. This model accounts for ~ 17 per cent of the FIRB. 



/fires = 5.0 x 1CT 10 Wm~ 2 sr _1 from Fixsen et al. (1998), 
we calculate the fraction of the FIRB which can be ac- 
counted for by our obscured QSOs. Low qo models provide 
~ 8 — 17 per cent of the FIRB intensity, compared with 
~ 3 — 8 per cent for the qo = 0.5 models, as shown in Tabic 

I 

The difference in the predicted number counts and back- 
ground intensity between the fcov = 1 and /« = 0.5 cases 
is due to the increased intrinsic far infra-red luminosity of 
sources which have an isotropic absorbing medium, as a 
greater fraction of the nuclear radiation is intercepted. These 
intrinsically brighter sources are therefore detected at higher 
fluxes, increasing the bright-end number counts. However, 
at faint fluxes, the predictions are similar for both assumed 
covering factors, implying that we are observing high red- 
shift objects, which by definition will be the most luminous 
sources in either case. 

For an isotropic covering medium, the predicted source 
counts of obscured AGN at 2 mJy account for between 8 and 
20 per cent of the total counts, whereas for the anisotropic 
case, the numbers drop by about a third, since the objects 
are less luminous. The contribution to the intensity of the 
FIRB is ~ 3 — 13 per cent in the f CO v ~ 0.5 case, which is 
about 25 per cent lower than the ~ 4— 17 per cent predicted 
for fcov = 1. 



4.3 The effects of changing z max and Td 

If we were to assume that obscured QSOs exist out to red- 
shift Zmax ~ 10, then we can make a much larger contri- 
bution to the source counts and FIRB intensity, as shown 
for Td = 30 K and f CO v — 1 in Fig. ^|(a). However, it is 
difficult to envisage a scenario in whic h large numbers of 
dusty QSOs form at such early epochs ( |Efstathiou fc Reel 
1985), and therefore this places a useful constraint on the 
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Figure 6. (a) The effect on the observed 850^tm source counts, 
of varying the maximum rcdshift, -Smaxi 

for dust at T d = 30 K, 
for z max = 2, 5 and 10. (b) The effect on the observed 850/im 
source counts, of varying the dust temperature between Td = 
30 K and T d = 70 K, over redshifts < z < 5. In both cases, 
parameters for our fiducial model POW(O.O) have been used, for 
seven populations with a flat column distribution, and a covering 
factor f cov = 1. 



maximum contribution such models can provide. Obscured 
QS Os from < z < 10 could contribute over 30 per cent of 
(1998) measurement of Jpirb- 



Fixsen et al 



the 

There are very few complete sub-mm surveys of AGN at 
low and high redshifts in the literature, with enough data to 
constrain the dust temperatures in such objects. The uncer- 
tainties involved in making such estimates are high, where 
assumptions have to be made about the cosmology, the dust 
masses and the opacities, and in addition, observations at 
several different wavelengths are required in order to start 
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Figure 7. The predicted total source counts for the qo = 0.5 
models, for power-law (solid lines) and polynomial (dashed lines) 
luminosity evolution, using a tilted distribution of column densi- 
ties. Predictions for both a covering factor of f cov = 1 (bold lines) 
and f cov = 0.5 (thin lines) have been plotted, using T d = 30 K 
and Zmax — — 5. 



to put limits on the parameters in the models. Treating the 
dust as isothermal is likely to be an oversimplification, but 
will be adequate for our purposes, and therefore we estimate 
that the range 30 K < Td < 70 K should span the most prob- 
able temperatures. 

For each column density used, the FIR luminosity is in 
effect fixed by the known amount of nuclear X-ray and op- 
tical luminosity which has been absorbed. If we increase the 
assumed dust temperature, the mass of the dust must be 
decreased in order to keep Lpir constant. As the tempera- 
ture is increased, the peak of the thermal radiation moves 
to higher frequencies, thereby reducing the intensity of the 
source at sub-mm energies. In Fig. ^jb), we show the effect 
of varying the temperature on the sub-mm counts, where 
the highest number of sub-mm sources are observed for the 
model with the lowest dust temperature. If the dust around 
most AGN is warm, Td ~ 70 K, then the contribution to 
the 850pm FIRB will be negligible, whereas at 100/im, their 
impact will be more significant. 



4.4 Tilted column distributions for qo = 0.5 models 

One further test here for the obscured QSO model is to 
look at the tilted distribution of column densities, which 
was invoked in order to obtain a better fit to the XRB for 
high density, qo = 0.5, models for the XLF. By skewing the 
distribution of objects towards higher column densities, this 
could have the effect of overpredicting the source counts and 
background intensity. However, in Fig. ^| it can be seen that 
the source counts predicted with this distribution for a range 
of models lie well below the observed counts, and flatten off 
towards fainter fluxes, again primarily due to volume effects. 
The maximum contribution to Jfirb from these models is 
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Figure 9. The predicted contribution to the FIRB from seven 
populations of obscured AGN with a flat column distribution, 
using power-law and polynomial luminosity evolution, for both 
qo = 0.0 and 50 = 0.5 cosmologies, as labelled. A covering factor 
fcov = 1 is used, with a dust temperature = 30 K, integrated 
over redshifts < z < 5. The data are as described in Section Ep 
CMB - Cosmic Microwave Background, Mather et al. (1994); F98 
- Fixsen et al. (1998); P96 - Puget et al. (1996); H98 - Hauser et 
al. (1998); S98 - Schlegel et al. (1998). 
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Figure 8. Predicted number-redshift distribution for an 850/^m 
survey to 2 mjy, of obscured QSOS with Tjj = 30 K, and z max = 
5, for covering factors of (a) f cov = 1, and (b) f cov = 0.5. Note 
the higher number densities predicted for low qo models. 



15 per cent, for model POW(0.5)t with T d = 30 K and 



fcov — 1 • 



4.5 Number-redshift distributions 

In Fig. ^, we plot the number-redshift distributions pre- 
dicted by our obscured QSO models, for an 850/zm survey 
to 2mJy (c.f. the la confusion limit for JCMT/SCUBA of 
0.44 mjy beam - ; Blain, Ivison & Smail 1998). Again, we 
use a dust temperature of Td = 30 K, and consider sources 



out to ledshift z ma x = 5, and make predictions for both (a) 
fco V = 1 and (b) f cov = 0.5. We have taken four models with 
a fiat column distribution of seven populations of obscured 
AGN, and two models using a til ted distribution of six pop- 

(all for f cov = 0.5). In 



ulations, as described in Section 3.1 



the go = 0.0 cosmology, the number of sub-mm sources pre- 
dicted at redshifts z > 3 is much higher than for go = 0.5, 
consistent with the steeper faint-end slope of the LogN:LogS 
relation for a low go universe. Once more source identifica- 
tions become available, it will be very straightforward to 
determine the form of the luminosity evolution, due to the 
fact that it is as easy to detect a sub-mm source at high 
redshift as at low redshift, and therefore if large numbers 
of high redshift obscured QSOs exist, they will be found in 
deep SCUBA surveys. 

The follow-up identification programs for the latest 



SCUBA sub-mm surveys are nearing completion, and the 
source catalogues will soon be published. We shall there- 
fore be able to compare the numbers of AGN detected with 
the redshift distributions predicted by our models. It is in- 
teresting to note that ~ 10 per cent of the optical coun- 
terparts to sub- mm sources are clas sified as Extremely Red 
Objects (EROs; Hu & Ridgway 1994 ) , wh ich are found in K - 
band images of the SCUBA error boxes dSmail et al. 1999[ ). 
Deep radio maps provide more accurate positional informa- 
tion than the SCUBA maps, assuming that the radio and 
sub-mm emission is due to the same mechanism, and the 
ERO counterparts are confirmed by the radio data. It will 
be very difficult to obtain spectroscopic redshift informa- 
tion about these objects, even with a 10-m class telescope, 
as they are so faint in the optical, with I S> 25. However, 
near infra-red spectroscopy has been used success fully to ob- 



tain a redshift of z = 1.44 for the ERO HR10 (Deyetal 



199£). Potentially, these could be examples of very highly 



obscured AGN at high redshift, for which in the optical and 
near infra-red we see the dusty host galaxy (hence the very 
red colours), and only see evidence for the AGN through 
the sub-mm and radio emission. Future deep surveys with 
AXAF and XMM, such as the proposed AXAF and SCUBA 
observations of the SSA13 field (Cowie et al.), will provide 
a vital test of these theories, in searching for faint X-ray 
emission associated with these sources. However, unlike in 
the sub-mm regime, the X-ray fc-correction is not working in 
our favour, and if we assume that these sources are at high 
redshift, then this is still an extremely ambitious project. 
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4.6 Predicted spectrum of the far 
infra- red/sub-mm background 

We have shown from the sub-mm number counts predicted 
by our models that obscured AGN provide a small but non- 
negligible fraction of the intensity of the far infra-red back- 
ground at 850/^m, ranging from 1 to 33 per cent. Fig. ^Jshows 
the spectrum predicted from our flat distribution of column 
densities using a covering factor f cov = 1, a dust tempera- 
ture Td = 30 K, and integrated over redshifts < z < 5. It 
can be seen that all four models predict the same intensity 
at high frequencies. This is to be expected, as the emission 
here is dominated by low redshift objects for which the rest- 
frame peak of the 30 K thermal spectrum is at ~ 100/im. 
Since the model parameters have been obtained from fits of 
the X-ray luminosity function to X-ray selected QSOs from 
ROSAT and the Einstein EMSS, which have relatively low 
median redshifts of z ~ 1.5 and z ~ 0.2 respectively, we 
would expect the predictions to diverge at high redshift and 
therefore low frequencies. 

In Fig. [h| we show the effects of changing the maxi- 
mum redshift and the dust temperature used in the models. 
In panel (a), it can be seen that the low redshift sources 
contribute at higher frequencies around 100/im, whereas the 
high redshift sources account for the largest fraction of the 
FIRB intensity at 850/im. In panel (b), we show how as 
the dust temperature is increased, the peak of the predicted 
FIRB background moves to higher energies. The integrated 
luminosity is fixed by the amount of absorbed nuclear X- 
ray and optical radiation, and is independent of tempera- 
ture. Since this is a spectral energy distribution diagram, in 
which equal areas mean equal energies, then the predicted 
intensities for each model are identical after a shift along 
the frequency axis. A dust temperature of Td ~ 40 K would 
give the maximum contribution to the peak of the observed 
FIRB. 



5 DISCUSSION 

The first question to address is what these obscured AGN 
will look like at wavelengths other than the sub-mm. Com- 
pared with the majority of the X-ray source population, sub- 
mm sources have much higher dust masses and therefore 
much higher column densities. Hence in general, the opti- 
cal and X-ray nuclear emission will be totally obscured, and 
the counterparts are likely to appear like relatively "normal" 
galaxies. These galaxies may perhaps have narrow emis- 
sion lines in their optical spectra, possibly of high ioniza- 
tion species, or may look very dusty from their optical and 
infra-red colours. However, if the obscuring material is not 
assumed to be isotropic, then a fraction of these highly ob- 
scured sources will be orientated such that our line of sight 
lies within the opening angle of the torus, and therefore 
the X-ray and optical nuclear flux will escape unattenuated, 
while simultaneously a large sub-mm flux is detected from 
the dust in the torus. 

The identification of optical counterparts to sub-mm 
sources has similar problems to those encountered with X- 
ray data, in that the point-spread function of the telescope 
is large (Half Power Beam Width = 14"7 at 850^m for 
SCUBA/JCMT), and therefore a number of plausible coun- 
terparts can lie in the error box. As sub-mm sources often 
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Figure 10. (a) The effect on the predicted FIRB spectrum of 
varying the maximum redshift, Zmaxi 

for dust at T d = 30 K, 
for z m ax = 2, 5 and 10. (b) The effect on the predicted FIRB 
spectrum of varying the dust temperature between T<j = 30 K 
and Td = 70 K, over redshifts < z < 5. In both cases, our 
fiducial model has been used, for seven populations with a flat 
column distribution, and a covering factor f cov = 1. The data are 
as described in Section ^| 



have associated radio emission, deep radio m aps of sub-mm 



survey fie lds have been taken (tvison et al. in preparation 



Richards 1998) . Star- formation regions are expected to con- 
tain many supernova remnants, known to be strong radio 
emitters, and even radio-quiet AGN are likely to be detected 
as very faint, /j,3y sources in extremely deep radio maps. The 
high angular resolution of the radio data combined with the 
fact that radio emission is not affected by the presence of 
dust, means that optical counterparts to the sub-mm sources 
can be found in a much less ambiguous manner. 

Once the counterpart has been found, the mechanism 
giving rise to the sub-mm emission must then be determined. 
By assuming that for starburst galaxies, the star-formation 
rate controls b oth the radio e mission and the thermal sub- 
mm emission, Carilli & Yun (1999) use the radio to sub- 



mm spectral index as a redshift indicator. If however, an 
independent measure of the source redshift can be obtained, 
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then the radio to sub-mm spectral index can be used to 
infer whether there is any additional contribution to either 
component due to the presence of an AGN. A radio-loud 
AGN will have proportionally higher radio emission than a 
starburst galaxy, whereas an obscured radio-quiet AGN will 
have lower radio and higher sub-mm emission. 

At present, there is little information available with 
which to constrain the dust temperature or range of tem- 
peratures that should be adopted for our models. Here, we 
have assumed that Td remains constant with redshift, but as 
shown in Fig. || the temperature affects the predicted source 
counts significantly. The advent of larger samples of AGN at 
both low and high redshifts with multi- wavelength FIR/sub- 
mm data will enable any evolution of the dust temperature 
to be constrained. 

Having chosen to keep Td constant with redshift, and 
a universal opacity law, the only unconstrained parameter 
defining the luminosity is therefore the dust mass, Ma, (see 
Eqn [l]). Therefore, since the FIR luminosity is directly re- 
lated to the X-ray luminosity, which we have modelled to 
undergo pure luminosity evolution with redshift, then one 
possible interpretation is that the dust mass also evolves 
with redshift, with the form Md oc (1 + z) 3 . If the dust mass 
does not scale with the luminosity, then in order for our 
models to hold, there must be a balance between dust mass 
and temperature, with the mass compensating for changes 
in temperature. 

A further consideration is the question of whether the 
assumption of constant gas to dust ratio at all redshifts is ap- 
propriate. For a QSO, there exists a radius within which the 
radiation field is so intense that dust particles will not sur- 
vive, the so-called dust sublimation radius. Granato, Danese 
& Franceschini (1997) proposed that a significant quantity 
of gas exists inside the dust sublimation radius, and that 
a large proportion of the photo-electric absorption occurs 
within this region. This would have the effect of reducing the 
dust masses calculated using the column densities inferred 
from the observed photo-electric absorption at X-ray ener- 
gies. This would in turn, lower the sub-mm fluxes from ob- 
scured AGN, cutting their contribution to the source counts 
and the intensity of the FIR background. However, the cal- 
culation of the g as to dust rat io for SMM 02399-0136 at red- 
shift z = 2.8 by |Frayer et al] (1998) and SMM 14011+0252 
at redshift z = 2.6 by |Fraycr et at| (1999), from measure- 
ments of CO line emission combined with the sub-mm flux, 
gives a value similar to that found in nearby galaxies, from 
which they infer that CO emitting sources at high-redshift 
have already undergone significant chemical evolution. We 
have therefore chosen to adopt a constant gas to dust ratio 
in our models, since at present the available observations are 
insufficient to constrain any suitable alternatives. 



6 FUTURE OBSERVATIONAL TESTS OF THE 
MODELS 

The next decade promises to bring enormous advances in 
this field, with the advent of innovative instrumentation 
combined with the light-gathering capacity of 10-m class 
telescopes, plus the new generation of satellite-bourne detec- 
tors. In the near infra-red, large area surveys will be possible 
with new wide-field cameras, such as the Cambridge Infra- 
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Figure 11. Predicted number-redshift distribution for an 850/^m 
survey to 8 mjy, of obscured QSOS with = 30 K. and z max = 
5, for (a) f cov = 1 and (b) f cov = 0.5. As in Fig. pi we plot the 
qo = models (solid lines) and qo = 0.5 models (dashed lines) 
for seven populations with a flat column distribution, and the 
tilted qo = 0.5 models with six populations (dotted lines). The 
power-law models are plotted with bold lines in each case, and 
the polynomial models are plotted with thin lines. 



Red Survey Instrument (CIRSI; Beckett et al. 1997), allow- 
ing follow-up observations of the entire survey area for exist- 
ing and future X-ray surveys, rather than the pointed obser- 
vations of individual sources used previousl y. In the far infra- 
red, the FIRBACK survey ( [Lagache 1998| ) is an ISOPHOT 
175pm survey of 4deg 2 of sky at high Galactic latitudes, 
with the aim of determining the source populations making 
up the far infra-red background. Proposed satellite-borne 
mid and far infra-red observatories in clude the NA SA Space 
Infra-Red Telescope Facility (SI RTF; |Werner 199§| ), and the 
ESA Far Infr a -Red and Sub -millimetre Telescope (FIRST, 
Pilbratt 1993[ |Genzel 1997j ) scheduled for launch in 2001 
and 2007 respectively. 

In the sub-mm regime, a proposed wide area survey 
with SCUBA plans to cover a subset of the European La rge 
Area ISO Survey region (ELAIS; pliver et al. 1998[ ) of 
around 640 square arcminutes (0.178 deg^). The survey will 
have a brighter flux limit than existing deep pencil-beam 
surveys, and aims for a 3a detection threshold of ~ 8 mjy, 
from which ~ 40 sources are expected. Our predictions for 
the number-redshift distribution for such a survey are pre- 
sented in Fig. llKa) for f cov — 1, with the more conservative 
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predictions using f cov = 0.5 in Fig. 0(b). The expected REFERENCES 



number of sources ranges between ~ 14deg" 2 (POL(0.5), 



fcov 



0.5, 



5) and ~ 180deg~^ (POW(0.5)t, 



1999, MNRAS, 305, L59 



fcov = 1, z m ax = 5), depending on the evolution and col- 
umn density distribution used in the model. Clearly, if the 
obscuring torus is also a site of active star-formation and the 



sub- mm flux of each component is comparable (Frayer et al 
199 8j), j hen the number of AGN may be somewhat higher 
than predicted on the basis of the above model. 



7 CONCLUSIONS 

In this paper, we have exte nded the obscured QSO model 
for the X-ray background of Sunn fc Shanks (1999) to the 



sub-millimetre regime, by considering the fate of the X-ray, 
ultra-violet and optical energy absorbed by the obscuring 
medium. This energy goes into heating up the dust in the ob- 
scuring material, which then radiates thermally at far infra- 
red and sub-millimetre wavelengths. We have modelled the 
obscuring medium as either isotropic or having a toroidal 
geometry, which then dictates the intrinsic column density 
distributions which are consistent with the line-of-sight dis- 
tributions found from X-ray and optical observations. 

Since the spectrum in the sub-mm is rising steeply to 
higher frequencies, the fc-correction obtained is such that our 
obscured QSOs are equally visible for redshifts of 1 < z < 
10. We therefore use the observed sub-mm source counts at 
850/im and the spectrum of the far infra-red background to 
constrain our models, by ensuring that the large quantities 
of cool dust invoked would not exceed the observed emission. 

We have shown that a variety of plausible obscured 
AGN models, which provide good fits to the X-ray back- 
ground spectrum and the number counts at soft and hard 
X-ray energies, are consistent with the observed sub-mm 
source counts and intensity of the FIRB. The models pre- 
dict between 1 and 33 per cent of the FIRB intensity, and a 
similar fraction of the number counts, depending on how ex- 
treme is the model, but with the more conservative models 
predicting between 5 and 15 per cent. This is in good agree- 
ment with the fact that the majority of sub-mm sources are 
identified as starburst galaxies, providing the complemen- 
tary sources. In addition, the obscured AGN models may be 
a suitable candidate for the source of the sub-mm emission 
associated with EROs. 

Finally, we have made predictions of the redshift dis- 
tribution of sub-mm sources from our models, for both ex- 
isting and proposed surveys, and described how combined 
X-ray and sub-mm survey data will be able to determine 
the extent of the obscured QSO contribution to cosmologi- 
cal backgrounds at both high and low energies. 
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